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HIGHLIGHTS 


. We test a thermally powered circulator based on two-phase thermofluidic oscillations. 

• Two heat exchanger designs and two power transmission lines are compared. 

• Increased fluid inertia in the load improves the maximum flow-rate and stalling head. 

• The pumping capacity depends on the heat uptake and the liquid oscillation amplitude. 
. The stalling head is correlated to the internal working fluid temperature amplitude. 
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A thermally powered circulator based on a two-phase thermofluidic oscillator was constructed and 
operated successfully as a replacement for a central heating hot water circulator coupled to a domestic 
gas-fired boiler. During regular operation the thermally powered circulator demonstrated a pumped 
flow-rate that decreased monotonically as the head applied across it increased. A maximum measured 
flow-rate of 850 L/h was achieved at zero head, and a maximum head of 8.4 mH 2 0 was attained at 
near-stalling (zero flow-rate) conditions. In agreement with previous modelling studies of the technol¬ 
ogy, increased inertia in the load line seems to lead to improved circulator performance. Further, the 
oscillating circulator exhibited an operational frequency between 0.24 and 0.33 Hz, which was mostly 
determined by the circulator configuration. The pumping capacity was strongly affected by the oscillating 
liquid amplitudes in the power cylinder that defined the positive displacement amplitudes of the liquid 
piston into and out of the hot water circuit. The best circulator configuration was associated with lower 
operation frequencies and relatively large ratios of suction to discharge displacement. 

Crown Copyright © 2013 Published by Elsevier Ltd. All rights reserved. 


1. Introduction 

The depletion of finite fossil-fuel resources together with con¬ 
cerns, both health-related and environmental, relating to the re¬ 
lease of product gases from fossil-fuel combustion (CO x , NO x , 
SO x , particulates, etc.) to the atmosphere are acting as major driv¬ 
ers towards a higher efficiency and more sustainable energy future 
[1], Essential components of such a future, depending on the en¬ 
ergy demand in a particular setting, are the recovery and re-use 
of waste and low-grade heat, or the conversion of this heat to use¬ 
ful power. Low-grade (i.e. low temperature) waste heat is abun¬ 
dantly available in the industrial sector, and also in the domestic 
environment, while non-concentrated solar thermal and low-tem¬ 
perature geothermal energy are additional untapped energy 
sources with great potential if they can be harnessed simply, reli¬ 
ably and at low cost. 


* Corresponding author. Tel.: +44 20 759 41601. 
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The desire to utilise renewable energy resources where possible 
has led to an enhanced interest in thermoacoustic heat engines. 
Importantly, these systems are capable of operating directly with 
thermal energy inputs, such as solar or waste heat, and further 
have no, or few, moving components in addition to typically 
employing an inert gaseous working fluid with little or no adverse 
environmental impacts. The first thermoacoustic devices that dem¬ 
onstrated the possibility of heat conversion to sound (acoustic 
power) were the Sondhauss and Rijke tubes [2], The detailed study 
of the thermoacoustic oscillations in these devices has formed a 
basis on which thermoacoustic heat engines have been developed, 
by amongst other, Wheatley et al. [3,4], Backhaus and Swift [5], 
Ceperley [6], Sun et al. [7], Qiu et al. [8], Redlich and Berchowitz 
[9], Dai et al. [10], Luo et al. [11] and Tang et al. [12]. Engines based 
on travelling waves (in which pressure and specific volume are 90° 
out of phase) with reversible heat transfer have been shown to out¬ 
perform standing wave based systems (in which pressure and spe¬ 
cific volume are in phase) that require intrinsically irreversible 
heat transfer to operate [5[. 


0306-2619/$ - see front matter Crown Copyright © 2013 Published by Elsevier Ltd. All rights reserved. 
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Nomenclature 





ATjnt 

difference between r HH x-fiuid and T CH x-nuid (°C or 1<) 

Abbreviation 

^boil-in 

temperature of cooling flow driven by the electric pump 

NIFTE 

Non-Inertive-Feedback Thermofluidic Engine 


entering the boiler (°C or K) 

CH 

Central Heating 

PH20 

density of water (kg/m 3 ) 

TPC 

Thermally Powered Circulator 

Cp,H20 

specific heat capacity of water at constant pressure (J/ 

HHX 

Hot Heat Exchanger 


kg K) 

CHX 

Cold Heat Exchanger 

y + 

upwards/suction stroke amplitude of oscillation in the 

PTL 

Power Transmission Line 


power cylinder (m) 

HP 

High Pressure 


downwards/power/discharge stroke amplitude of oscil- 

LP 

Low Pressure 


lation in the power cylinder (m) 

TC-# 

Thermocouple 

Ap 

cross-sectional area of the power cylinder (m 2 ) 

T-# 

Temperature Dials 

f 

operation frequency of the TPC (Hz) 

HPR 

High Pressure Reservoir 

fo 

maximum f at zero head conditions (Hz) 

LPR 

Low Pressure Reservoir 

po 

density of water at 20 °C (kg/m 3 ) 

DAQ 

Data Acquisition 

Ps 

density of stainless steel at the actual measuring tern- 

SC 

Signal Conditioning 


perature (kg/m 3 ) 

PG-# 

Pressure Gauge 

p 

density of water at the actual measuring temperature 

rms 

Root Mean Square 


(kg/m 3 ) 



t 

Time (s) 

Symbols 




U 

flow-rate in the TPC hydraulic circuit (L/h) 

Subscript 


U e 

flow-rate in the CH hydraulic circuit (L/h) 

0 

zero head conditions; or 20 °C conditions (for rotameter 

U 0 

maximum flow-rate in the TPC circuit at zero head con- 


correction) 


ditions (L/h) 

s 

stainless steel /solid (for rotameter float material) 

A P 

pressure difference across the TPC (mH 2 0) 

max 

stalling, or near-stalling conditions 

A P 

time-averaged AP across the TPC (mH 2 0) 

p 

power cylinder 


maximum AP across the TPC at near-stalling conditions 

boil 

boiler 


(mH 2 0) 

H.loss 

heat losses from the boiler faceplate 

Qboil 

heat from the flame generated by burning natural gas in 

H 2 0 

water 


the boiler (W) 

int 

internal to the TPC 


heat losses from the boiler faceplate based on direct 

e 

central heating circuit, or external 


measurements (W) 

HHX-fluid fluid in the HHX 

Qhhx 

heat supplied to the HHX (W) 

CHX-fluid fluid in the CHX 

Qch 

heat supplied to the CH circuit (W) 

boil-in 

cooling flow driven by the electric pump entering the 

i-PTL 

length of the PTL (m) 


boiler primary CH heat exchanger 

^HHX-nuid 

temperature of the steam at the top of the power cylin- 

boil-out 

cooling flow driven by the electric pump exiting the boi- 


der at the moment of evaporation within the HHX (°C or 


ler primary CH heat exchanger 


K) 

P 

pressure 

TcHX-fluid 

temperature of the liquid water inside the CHX at the 

CH 

central heating, or central heating circuit 


moment of condensation (°C or K) 



^boil-out 

temperature of the cooling flow driven by the electric 




pump exiting the boiler (°C or I<) 




Fluidyne engines are a type of travelling wave thermoacoustic 
engine that employ liquid pistons thereby eliminating the need 
for mechanical moving parts and seals. A significant effort has gone 
into their development for fluid pumping applications [13,14]. The 
efficiency of a Fluidyne engine decreases at high loads, i.e. while 
pumping against large pressure differentials, due to high frictional 
losses in the tuning line, which are required for the satisfactory 
operation of the engine due to its inherent reliance on inertia 
(i.e. fluid mass). By partially evaporating the liquid displacer piston 
[15,16], higher power densities and larger displacement ampli¬ 
tudes at very low temperature differences were achieved, while 
paying a penalty in the form of poor thermal efficiencies. 

This low thermal efficiency at high loads was attributed to the 
excessive power required to drive the fluid oscillations in the tun¬ 
ing line, which motivated efforts to devise alternate systems that 
eliminate the dependence on inertia [17], This view has led to 
the introduction of a technology referred to as the ‘Non-Inertive- 
Feedback Thermofluidic Engine’ (NIFTE). This promising technol¬ 
ogy was recently proposed in Smith [17], who built an early 
small-scale prototype that employed an electrical heater as a heat 


source and devised a simple dynamic model for the purposes of 
early-stage engineering design. The NIFTE has been shown to be 
well suited to the conversion of low-grade heat to hydraulic (fluid 
pumping) work, and has been demonstrated to operate across tem¬ 
perature differences between the heat source and the ambient of as 
low as 30 °C [18]. 

The NIFTE is as a two-phase thermofluidic oscillator, an oscilla¬ 
tory heat engine in which persistent thermodynamic oscillations 
(e.g. of pressure, volume and temperature) are generated and sus¬ 
tained within a connected space of chambers and tubes by steady 
external temperature differences. The oscillations are driven by 
and give rise to heat and fluid flows, which involve the periodic 
evaporation and condensation of the working fluid, which is con¬ 
tained within the device simultaneously in both phases. 

With its few moving parts compared to mechanical heat en¬ 
gines as well as its lack of reliance on lubricants and dynamic seals, 
in common with most gas-phase thermoacoustic heat engines, the 
NIFTE promises improved reliability, extended lifetimes and low 
maintenance costs, while at the same time high manufacturing tol¬ 
erances and a lack mechanism may (potentially) allow reduced 
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capital costs. On account of these characteristics the NIFTE has at¬ 
tracted attention, specifically with regards to the performance 
capabilities of this technology. 

The NIFTE has been the subject of a number of follow-on studies 
[18-22] whose aim has been to try to gain a better understanding 
of the coupling between the underlying thermal, fluid flow and 
thermodynamic processes and how this gives rise to the dynamic 
behaviour (i.e. oscillatory operation) of the device. In these studies, 
a number of first-order models for the NIFTE were proposed with 
emphasis on the improved description of key components and pro¬ 
cesses in the device. However, there has been a lack of experimen¬ 
tal data against which these models can be validated. In this paper 
we provide experimental results of the NIFTE as a replacement for 
a central heating hot water circulator in a domestic setting. It is 
significant to point out that the present paper constitutes the first 
report on the employment of an early prototype of this technology 
in a practical application. 

2. Experimental methods 

2.1. Apparatus 

A schematic representation of the experimental setup em¬ 
ployed in the present study is shown in Fig. 1. The apparatus was 
designed so as to allow the flexible variation of the input variables 
(i.e. the fuel burn rate in a domestic boiler burner, the flow-rates of 
fluid through the flow circuits), as well as access for temperature, 
pressure and flow-rate measurements at key positions in the sys¬ 
tem. Two separate fluid circuits were involved in the arrangement: 
(i) the central heating (CH) circuit, described in Section 2.1.1 below 
and (ii) the NIFTE [17-22] thermally powered circulator (TPC), de¬ 
scribed in Section 2.1.2. Water was chosen as the preferred liquid 
for both flow circuits. In the present experimental investigation 


the two circuits were kept separate in order to allow a more de¬ 
tailed understanding of the overall system by imposing indepen¬ 
dent variations to the two circuits, for example by independently 
varying the two flow-rates and pressures. It is intended that in a 
later stage the two circuits will be combined, thus creating a single 
central heating loop driven by the TPC. 

2.1.1. Central heating (CH) system and boiler 

The TPC was attached to a modern, commercially available cen¬ 
tral heating system (Vaillant ecoTEC 824), modified for the pur¬ 
poses of the current investigation. The ecoTEC 824 is a wall-hung 
domestic ‘high-efficiency’ (~90%) condensing combination boiler. 
Condensing boilers recover waste heat from the flue gases for in¬ 
creased efficiency, while combination boilers are capable of both 
central heating and direct hot water supply. 

A pre-programmed fan unit within the boiler mixed natural gas 
with laboratory air. The resulting fuel-air mixture was sent to the 
combustion chamber where it was burnt in a spark-ignited pre¬ 
mixed flame. Fig. 1 shows that the flame was stabilised on a perfo¬ 
rated cylindrical surface (0). The heat from the flame Qboii was used 
to heat the water (1) in the central heating (CH) circuit directly, as 
well as to provide heat to the evaporator or ‘hot heat exchanger’ 
(HHX; 9) of the TPC (this is described in detail below, in Sec¬ 
tion 2.1.2). The boiler heat setting was adjusted by an automatic 
control valve that regulated the flow-rate of the fuel-air mixture 
burnt in the combustion chamber, generating heat between 7 
and 19 kW. 

The hot (40-70 °C) CH water was pumped out of the boiler by 
an electric circulator pump (5) supplied with the boiler. Typically 
in a domestic central heating setting, the CH water is then 
circulated for space heating through a number of radiators 
placed throughout the household. In the present experimental 
arrangement the CH water was circulated instead, first through a 



Fig. 1. Schematic of the apparatus and measurement facility. Showing: (0) the flame stabilised on the perforated cylindrical surface of the burner in the boiler; (1) the water 
flow in the central heating (CH) circuit; (2) the rotameter for the measurement of the CH flow-rate; (3) the CH circuit expansion vessel; (4) the large fan-assisted heat 
exchanger unit for heat removal; (5) the CH circuit electric circulator pump; (6 and 7) the inlet and outlet of the CH flow through the outside (shell-side) flow channel of the 
cold heat exchanger (CHX); (8) the condenser or CHX; (9) the evaporator or hot heat exchanger (HHX); (10) the power cylinder of the NIFTE thermally powered circulator 
(TPC); (11) the ‘HHX-CHX connection’ between the HHX and the inside (tube-side) channel of the CHX; (12) the ‘vapour connection’ between the HHX and the TPC power 
cylinder; (13) the ‘power transmission line’ (PTL) that connected the bottom of the power cylinder to the TPC flow circuit; (14) the ‘feedback line’ between the TPC power 
cylinder and the CHX; (15) the Y-fitting in the flow circuit of the TPC; (16 and 17) the two check (one-way) valves in opposite orientations in the TPC flow circuit; (18 and 19) 
the two vertical cylindrical reservoirs/pressure accumulators in the TPC flow circuit; and, (20) the combination of a ball and gate valve placed in series to control the flow-rate 
and pressure in the TPC flow circuit. Also showing the flow-rate measurement in the CH circuit U e and the flame in the boiler burner (with heat setting Qboii) from which heat 
Qch is supplied to the CH circuit and heat Qhhx to the HHX of the TPC. Temperatures TC-1 and TC-2 were sampled continuously, while T-l and T-2 were recorded from 
displays at steady state. 









































C.N. Markides, A Gupta/Applied Energy 110 (2013) 132-146 


135 


rotameter (2) and an expansion vessel (3), and then through a large 
fan-assisted heat exchanger unit (4) that acted as the heat sink, 
rejecting the heat into the laboratory and thus cooling the CH 
water flow. The laboratory was ventilated continuously with a 
large (300 W; 0.65 m 3 /s) extractor fan, thus maintaining the labo¬ 
ratory temperature at 29 ± 3 °C. The rotameter was used to mea¬ 
sure the flow-rate in the CH circuit. The indicated flow-rate read 
from the scale was corrected for density variations due to temper¬ 
ature using the standard multiplier correction term [p 0 (Ps - p)l 
P(Ps - Po)] 1/2 l where p 0 = 998 kg/m 3 and p are the water densities 
at the calibrated (20 °C) and the actual measuring conditions 
respectively, while p s = 8030 kg/m 3 is the density of the stainless 
steel float. The maximum applied correction was 1.2% at 70 °C. 
The expansion vessel was used as an accumulator, to steady the 
pressure and flow in the CH circuit. Specifically, it was used to ab¬ 
sorb any pressurisation due to thermal expansion of the liquid 
water within the relatively short CH pipe work, but also in the case 
of steam formation in the boiler at low CH water flow-rates. 

After the fan-assisted heat exchanger unit, the cold (30-40 °C) 
water was pumped back into the boiler by the electric pump, via 
the outside (shell-side) cooling flow channel (6 and 7) of the con¬ 
denser or ‘cold heat exchanger’ (CHX; 8). In Fig. 1, (6) denotes the 
inlet and (7) the outlet of this stream. Thus, the flow-rate in the 
closed-loop CH circuit U e was also affected by the flow resistance 
in the CHX of the TPC, and specifically in the shell-side cooling flow 
channel as shown in Fig. 1. Two different CHX designs were used 
(described in Section 2.1.2), which resulted in CH water flow-rates 
of U e = 250 and 320 L/h depending on the configuration (refer to 
Section 2.2 for further details). The electric pump was always run 
at full power. 

In summary, the TPC being essentially a heat engine, rejected its 
heat initially to the CH water loop at the CHX, resulting in some 
preheating of the CH water before this entered the boiler, after 
which it was rejected to the laboratory by the fan-assisted heat ex¬ 
changer unit, and finally to the atmosphere by the extractor fan. 
The level of pre-heating was measured at between 5% and 35% of 
the combustion chamber generated Qboii in the boiler. 

2.3.2. Thermally powered circulator (TPC) 

Referring again to Fig. 1, the thermally powered circulator (TPC) 
comprised three main components: (i) the hot heat exchanger 
(HHX; 9); (ii) the cold heat exchanger (CHX; 8); and (iii) the power 
cylinder (10). The HHX (9) consisted of tubes introduced into the 
boiler combustion chamber. At their lower end (11) the HHX tubes 
were connected to the inside channel (tube-side) of the CHX and at 
the top (12) they were connected to the top of the power cylinder. 
We will refer to these as the ‘HHX-CHX connection’ and 'vapour 
connection’, respectively. 

All HHXs and CHXs were custom made. Two different HHX de¬ 
signs were used to match two different designs for the CHX: (i) a 
large-capacity heavily finned HHX with large diameter channels 
(‘TPC Configuration A’) and (ii) a less intrusive HHX design with a 
reduced fin area and smaller diameter channels (‘TPC Configura¬ 
tion B’). The corresponding designs for the matching CHXs were: 
(i) a long single shell-single tube concentric tubular CHX and (ii) 
a short single shell-multiple tube shell-and-tube CHX. The long 
tubular CHX comprised a single large diameter tube, whereas the 
short shell-and-tube CHX comprised a number of small diameter 
tubes in parallel. We will refer to the former as the ‘large-capacity’ 
HHX-CHX combination or TPC configuration, and to the latter as 
the ‘small-capacity’ HHX-CHX combination or TPC configuration. 

The power cylinder (10) was a circular tube constructed from a 
clear high-temperature plastic, which permitted the visual obser¬ 
vation of the internal vertical oscillations of the liquid. The choice 
of material also provided superior insulating behaviour compared 
to glass [23,24]. The cylinder was connected at the top to the 


HHX tubes via the ‘vapour connection’, while at the bottom it 
was connected to a liquid ‘power transmission line’ (PTL; 13) as 
well as to the tube-side channel(s) of the CHX. We will refer to this 
pump-side connection between the power cylinder and the CHX as 
the ‘feedback line’ (14). The feedback line between the CHX and the 
power cylinder was kept as short as possible, as allowed by the 
configuration and was not varied. 

The PTL connected the bottom of the power cylinder to the flow 
circuit (15) of the TPC. Two PTL hoses were used with different 
lengths Lpx L = 2 and 5 m (see Fig. 1), but which were otherwise 
identical. Experiments with the two lengths were done in both 
the large- and small-capacity configurations of the TPC. The TPC 
flow circuit included: a Y-fitting with two check (or one-way) 
valves (16 and 17), one on each branch but in opposite orienta¬ 
tions; two vertical cylindrical reservoirs (18 and 19), doubling as 
pressure accumulators (or expansion vessels); and, a combination 
of two valves placed in series (20), one ball valve and one gate 
valve. 

The flow generated by the TPC in the PTL was purely oscillating 
(zero-mean) [17-22]. This was transformed to a unidirectional 
flow with the use of the Y-fitting and check valves. When the oscil¬ 
lating TPC generated an outward displacement stroke (y~ > 0 in 
Fig. 1 ) and hence a discharge of liquid from the power cylinder into 
the power transmission line, the check valve on one arm of the Y- 
fitting (16) opened, sending liquid into the discharge or high-pres¬ 
sure (HP) reservoir (18). During this discharge process the other 
(reversed) check valve in the second arm of the Y-fitting (17) en¬ 
sured that liquid was not able to flow from the supply or low-pres¬ 
sure (LP) reservoir (19) through the Y-fitting back into the power 
transmission line. When the TPC generated an inward displace¬ 
ment stroke (y + > 0 in Fig. 1 ) and hence suction of liquid from the 
power transmission line back into the power cylinder, the check 
valve in first arm of the Y-fitting closed, and the check valve in 
the second arm of the Y-fitting opened allowing liquid to flow from 
the LP reservoir through the Y-fitting back into the power trans¬ 
mission line. The liquid that accumulated inside the HP reservoir 
was then forced through the ball and gate valves into the LP reser¬ 
voir, and finally through the second arm of the Y-fitting to close the 
TPC flow circuit (15). Thus, the pressure (or head) rise across the 
TPC was equal to the difference between the pressures in the HP 
and LP reservoirs. For the purposes of the present investigation, 
the ball and gate valves were used to adjust externally the pressure 
drop in the TPC flow circuit, and thus also the pressure difference 
across the TPC. 

Referring to Fig. 2, each reservoir was a vertically oriented cylin¬ 
drical volume filled partially with liquid and partially with air. The 
two connecting valves were adjusted to provide coarse and fine 
control of the mean pressure difference between the two reservoirs 
in the circuit, which was also the pressure rise across the TPC AP, 
and consequently also of the flow-rate from the HP to the LP reser¬ 
voir, i.e. the flow-rate in TPC flow circuit U. The LP reservoir was 


Pressurisation point 
Closed-Ji Open-^n 


PG-1 

TPC flow ( U) 
from 
check valve 1161 




^Scale 

HP 

— 

LP 

Res ' 181 


Res ' 191 

--- 

— 

-- 


Ball/gate 


> TPC flow ( U) 
to 


valves 


( 20 ) 


check valve 


( 17 ) 


Fig. 2. Pressure reservoirs and measurement details. 
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kept open to atmosphere, such that the pressure within it fluctu¬ 
ated slightly about atmospheric pressure. The pressure in the HP 
reservoir fluctuated about an elevated mean pressure determined 
by the valves, and was adjusted further by the addition or removal 
of compressed laboratory air from the top. The fluctuations in both 
pressures arose as a consequence of the pulsating displacement ac¬ 
tion of the TPC. The magnitude of these fluctuations was mini¬ 
mised effectively by maximising the percentage volume of the 
reservoirs occupied by air, resulting in amplitudes of no more than 
a few cm. 

Finally, the primary requirement from the working fluid of the 
TPC was that the chosen fluid evaporated in the HHX and con¬ 
densed in the CHX, subject to the temperatures within these com¬ 
ponents that could be achieved practically in our setup. A number 
of secondary requirements were also considered, such as cost and 
toxicity. The decision was made to select water as the working 
fluid for the TPC. Ultimately, the current study investigates the 
possibility of replacing the electric circulator in the CH circuit with 
the TPC, even though this has not been done here for measurement 
flexibility. An additional advantage of the decision to pick water as 
a common fluid in both the CH and TPC flow circuits is that any fu¬ 
ture attempt to connect the two circuits, thus removing the need 
for the electric circulator, will not require a device or scheme to 
maintain two different liquids separated. 

2.2. Measurement methodology 

Prior to operation, a note of the ambient temperature, pressure 
and humidity in the laboratory was made. A desired heat setting 
(fuel burn) Qboii was prescribed to the boiler control system, after 
which the CH boiler was switched on. Heat settings of Qboii = 7, 
10, 12 and 15 kW were used in the experiments presented in this 
paper. 

After switching on the boiler the fuel-air fan began to introduce 
the flammable mixture into the combustion chamber, while the 
electrical circulator began to circulate the water in the CH circuit. 
The water flow-rate in the CH circuit was U e = 250 L/h in all tests 
involving the large capacity TPC configuration (i.e. with the 
large-capacity HHX-CHX combination), and U e = 320 L/h in all 
tests involving the small-capacity TPC configuration (i.e. with the 
small-capacity HHX-CHX combination). A few seconds later con¬ 
tinuous sparking was generated in the combustion chamber until 
the flammable fuel-air mixture was ignited and a stable flame suc¬ 
cessfully established. As the TPC started up and the flow-rate U in 
the TPC flow circuit (and pressure in the HP reservoir) increased, 
the two coarse and fine flow control valves were adjusted to 
achieve a certain pressure difference AP across the TPC. It should 
be noted that this arrangement imposed a dynamic load on the 
TPC, that is, a pressure drop that resulted by a loss in dynamic 


pressure (or, dynamic head) in the flow through the valves in the 
CH flow circuit. This arrangement is different to that of a static load 
that would have been imposed had the TPC been connected so as to 
pump a liquid over a certain vertical height from one reservoir to 
another. 

Two types of measurements were made: (i) visual observations 
of certain quantities on dials, gauges or displays and recording 
after some time was allowed for the system (CH and TPC flow cir¬ 
cuits) to reach steady state and (ii) continuous measurements of 
important quantities with electronic instruments and sampling 
by a Data Acquisition (DAQ) system, from the moment sparks ap¬ 
peared in the boiler for fuel ignition. For the latter, a 16-bit Na¬ 
tional Instruments USB-6259 DAQ device was used to record the 
output from important thermocouples and all pressure sensors. 
Sampling was done at a frequency of 10 Hz. The raw sensor out¬ 
puts were sent to a National Instruments SC-2345 signal condition¬ 
ing connector block for low pass noise filtering, after which they 
where data-logged by the DAQ device onto a dedicated PC for pro¬ 
cessing. The DAQ system was switched on as soon as sparks ap¬ 
peared in the boiler combustion chamber, after which data was 
recorded over a time of 5 min. 

From this point, a set of variables was measured for each TPC 
configuration (i.e. choice of HHX, CHX and power transmission 
length Lptl). boiler heat setting Qboii, and control valve setting in 
the TPC flow circuit that effectively set the relationship between 
the flow-rate U and pressure difference AP. These included mea¬ 
surements at key positions: (i) within the TPC (the temperatures 
within the HHX, CHX, power cylinder, HHX-CHX connection and 
vapour connection): (ii) in the TPC flow circuit (the temperature 
and pressure inside the LP and HP reservoirs, and the liquid dis¬ 
placement amplitude and frequency in the power cylinder); and 
(iii) in the CH circuit (the circuit flow-rate U e , the temperature of 
the inlet and outlet flow from the CHX, and the temperature of 
the inlet and outlet flow from the boiler). These measurements 
are either reported directly, used to estimate dependent quantities 
that are presented in this paper, or used during the runs to confirm 
steady-state behaviour of the overall system. 

2.2.3. Pressure 

The pressure in the HP reservoir was measured with an elec¬ 
tronic pressure sensor (PG-1 in Fig. 2), the output from which 
was calibrated using a mercury manometer. The pressure in the 
LP reservoir was found to vary by no more than ±0.1 m about 
atmospheric pressure. The final pressures were further corrected 
for small hydrostatic differences in the reservoirs by using a scale, 
as indicated in Fig. 2. The reported mean pressure difference be¬ 
tween the HP and LP reservoirs, which was also the dynamic load 
pressure rise across the TPC AP, was averaged over 10 pulsations. 
For a typical measurement of AP see Fig. 3. In cases where the 




(a) (b) 


Fig. 3. Pressure difference measured across the TPC during (a) low-load, normal steady-state operation and (b) stalling due to pressure at high load. 
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Table 1 

Qhhx in kW at the various Qt, oil for different heat exchanger/thermally powered 
circulator configurations. 


Qhhx (kW) 


Qboii (kW) 




HHX-CHX/TPC 

configuration 

Lptl (m) 

7 

10 

12 

15 

Small-capacity 

2 

- 

2.1-2.9 

- 

4.4-4.7 


5 

- 

2.5-3.2 

- 

4.0-4.7 

Large-capacity 

2 

2.3-3.5 

3.0-3.6 

3.2-5.2 

- 


5 

- 

2.7-3.9 

4.1-4.7 

4.0-5.5 


pump stalled, the maximum attainable AP is reported. The total 
uncertainty in the reporting of AP is ±0.1 m. Also, it was estab¬ 
lished that cycle-to-cycle variations in AP had a standard deviation 
of 0.3 m. 

2.2.2. Temperature 

Two continuous temperature measurements were made using 
miniature (0.25 mm and 0.50 mm diameter) mineral insulated K- 
type thermocouples. The two thermocouples were placed: (i) at 
the top of the power cylinder (indicated by label ‘TC-1’ in Fig. 1) 
and (ii) inside the CHX near the HHX-CHX connection (‘TC-2’ in 
Fig. 1 ). These were used to report the mean temperature of the steam 
at the top of the power cylinder at the moment of evaporation within 
the HHX r HH x-fiuid. and the mean temperature of the liquid water in¬ 
side the CHX at the moment of condensation T C Hx-fiuid. respectively. 
From these two temperatures it is then possible to calculate the 
mean ‘internal’ temperature difference as seen by the thermody¬ 
namic cycle undergone by the working fluid of the TPC, i.e. AT in t = 
^HHx-fiuid - f'cHX-fiuid- The raw thermocouple signals were connected 
to an Audon TCK-4 amplifier, which also preformed on-board cold 
junction compensation. The output from the thermocouple ampli¬ 
fier was read by the DAQ system. Finally, the saved voltage data 
was converted to temperature based on information provided in 
the manufacturer’s datasheet and standard thermocouple conver¬ 
sion equations. 1 It has been established that T H Hx-fiuid has a total 
uncertainly of ±0.5 K and a cycle-to-cycle standard deviation of 
1.5 K, whereas T C Hx-nuid has a total uncertainly of ±1.0 K and a cycle- 
to-cycle standard deviation of 2.5 K. The total uncertainty includes 
calibration and non-linearity in amplification errors. 

Two further temperature measurements were made, indicated 
by labels ‘T-l’ and ‘T-2’ in Fig. 1. These were of the temperature 
of the CH flow entering Tb 0 n-in and exiting r b0 ii-out the boiler at stea¬ 
dy state. Both temperatures were read from an Omega HH-21 ther¬ 
mocouple display unit (stated accuracy from the manufacturer 
±1 K), such that the total error associated with these measure¬ 
ments is ±2 K. From these temperatures and the knowledge of 
the flow-rate (J e in the CH flow circuit it was possible to calculate 
the heat transferred from the flame into the CH circuit, Qc H = Ph 2 o 
Ue Cp,H 20 (Tboii-out - T b0 ii-in)- Further, given our knowledge of the to¬ 
tal flame heat Qboii, it was possible to estimate the heat from the 
flame that went into the TPC at the HHX, Qhhx = Qboii - Qch - 
Oh, loss, where Qn.ioss is an estimation of the heat losses from the 
boiler faceplate based on actual measurements of the faceplate 
temperature and area. The heat loss Qh.ioss was between 1.0 and 
1.7 kW, which was 10-15% Qboii, such that the ratio Qh.ioss/Qhhx 
was between 20% and 50%. It was higher at higher boiler settings 
Qboii due to higher faceplate temperatures. From Table 1 Qhhx 
was found to be between 2.1 and 5.5 kW. 

2.2.3. Flow-rate 

The vertical oscillatory motion of the liquid level (liquid-vapour 
interface) in the transparent power cylinder was observed by eye. 


1 Following http://srdata.nist.gov/its90/main. 


Upwards positive (suction)y + and downwards negative (discharge) 
y~ oscillation amplitudes were defined relative to an ‘equilibrium’ 
level that coincided with the mid-height of the HHX-CHX connec¬ 
tion, as illustrated in Fig. 1. The two amplitudes were read with the 
help of a ruler fixed directly onto the power cylinder in order to 
minimise errors. Multiple realisations were recorded for both y + 
and y~. In this paper we report amplitudes averaged over 10 real¬ 
isations, which were taken during a time when observations indi¬ 
cated that the TPC had settled to a steady state. 

For the frequency of oscillations we inspected the recorded 
pressure plots generated from the pressure gauge in the HP reser¬ 
voir (PG-1 in Fig. 2). We report an averaged frequency/calculated 
from the average period of oscillations over 10 steady-state cycles. 
The total uncertainly in the reporting of/is ± 0.01 Hz. The cycle-to- 
cycle standard deviation was 0.03 Hz. Hence, the TPC flow-rate was 
evaluated from U = (y + - y~)Apf, where A p is the cross-section of the 
power cylinder. 

Repeated runs were conducted in order to confirm that the per¬ 
formance of the TPC and the associated results were repeatable, 
but also to give an estimate of any variability. This process, to¬ 
gether with an error analysis based on each variable, revealed that 
the total uncertainly in U is ±15 L/h in the case of the small-capac¬ 
ity TPC and ±20 L/h for the large capacity configuration, respec¬ 
tively. However, as with other variables, cycle-to-cycle variations 
were more significant with a standard deviation of 50 L/h. 


3. Operational observations 

3.1. Normal operation 

The TPC is an oscillatory heat engine, whose operation relies on 
the continuous sequential alternation between two stages: (i) the 
power (or discharge) stroke and (ii) the suction stroke. In these 
stages positive displacement of fluid (in our case, water) is caused 
by the generation (by evaporation during the power stroke) and 
reduction (by condensation during the suction stroke) of a volume 
of water as this undergoes a two-phase thermodynamic cycle. A 
detailed description of the TPC cycle is available in any of the Refs. 
[17,21,22], 

Initially, with the boiler off, the entire internal volume of the 
TPC (including the HHX, CHX tubes, power cylinder and all connec¬ 
tions) contains water in the liquid phase at ambient temperature. 
Soon after the establishment of the flame in the boiler combustion 
chamber the heat originating from the hot combustion gases 
causes the water in the HHX to evaporate. After some initial tran¬ 
sients the TPC reaches steady-state oscillations. 

During the power stroke the HHX utilises a part (Qhhx) of the 
total heat (0^) generated in the boiler combustion chamber to 
evaporate water contained in the HHX. This occurs during the time 
when the water level (water-steam interface) moves vertically up¬ 
wards from the CHX into the HHX. The generated steam flows up¬ 
wards within the HHX tubes, through the vapour connection and 
into the top of the power cylinder. This causes the pressure in 
the power cylinder to increase and displaces the water in the 
power cylinder downwards and out into the HP reservoir. The 
effectiveness of the evaporation process determines the negative 
displacement power amplitude y~ in the TPC. The heat removal 
from the HHX due to the contact with cold water is associated with 
a quenching of the solid walls of the HHX tubes, leading to lower 
wall temperatures and progressively reduced steam generation 
ability in the HHX. Further, the steam introduced into the top re¬ 
gion of the power cylinder comes into contact with the colder walls 
in this space, which leads to partial condensation there. This loss 
mechanism, known as entrance condensation, was also identified 
in the early N1FTE prototype in Refs. [17,21], 
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(a) (b) 

Fig. 4. Temperature on the (a) HHX walls and (b) CHX walls, during normal steady-state operation. 


The power stroke process described in the above paragraph re¬ 
sults in a flow into and pressurisation of the HP reservoir. A clearly 
defined steam-water interface can be observed moving down¬ 
wards along the height of the power cylinder. The downward dis¬ 
placement of the water level in the power cylinder eventually 
results in a situation in which the water level drops below the 
water-steam interface in the HHX that is evaporating. This causes 
a hydrostatic pressure difference across the feedback line that con¬ 
nects the HHX with the power cylinder, which in turn eventually 
causes the water level in the HHX to drop into the CHX. 

The next stage is the suction stroke, which relies on the capacity 
of the CHX to remove heat from the steam contained within the 
CHX tubes. When the water-steam interface moves from the 
HHX into the CHX the steam condenses into water by contact with 
the cold walls. This leads to a reduced pressure (suction) that 
pulls the water in the power cylinder upwards and results in flow 
into the power cylinder from the LP reservoir. The steam-water 
interface can then be observed moving upwards along the height 
of the power cylinder. The effectiveness of the condensation stage 
determines the positive displacement suction amplitude y* of the 
TPC. When the water level in the power cylinder rises above the 
water-steam interface in the CHX that is involved in condensation 
there is a hydrostatic pressure difference across the feedback line 
that connects the CHX with the power cylinder. This pressure 
difference increases until it eventually causes the water level in 
the CHX to rise into the HHX, where the evaporation process 
associated with the power stroke will repeat itself. 

The heat rejected to the walls of the CHX during condensation 
in the suction phase flows into the CHX cooling water stream, 
and finally into the fan-assisted heat exchanger unit where it is 



(a) 


rejected to the laboratory. The heat deposited at the CHX increases 
the temperature of the CHX walls and progressively reduces its 
capacity for condensation. The durations of the power and suction 
strokes determine the oscillation period of the pump. 

Fig. 3a shows a typical time plot of the pressure difference 
across the TPC AP during normal operating conditions. This partic¬ 
ular low-load operating condition involves a near-zero time-aver¬ 
aged pressure difference AP. The oscillations are periodic, smooth 
and described relatively well by a single sinusoid, although the 
pointed peaks and troughs suggest that it does exhibit some (lim¬ 
ited) higher frequency content. This is an important point as it im¬ 
plies that, at least in these conditions, the TPC behaviour is 
dominated by processes whose frequency spectrum is discrete 
and characterised strongly by a single frequency component with 
some lower harmonics. 

Normal operation was associated with a root mean square (rms) 
in the oscillations of AP in the range 0.1-0.3 m (e.g. see Fig. 3a) for 
operation below AP = 1 m, and in the range 0.1-0.5 m for opera¬ 
tion up to AP = 5 m, increasing with increasing AP. This corre¬ 
sponded to a relative oscillation in the rms of AP with respect to 
the mean AP of up to 30% for AP < 1 m, to less than 20% for AP be¬ 
tween 1 and 6 m, to less than 10% for AP > 6 m, decreasing with 
increasing AP. 

The time plot in Fig. 3b was generated at high load and shows 
the behaviour of the TPC in stalling conditions. Note the three stall¬ 
ing and re-starting processes that have been captured during this 
recording. The stall points (at t ~25,120 and 265 s) are associated 
with high pressure and almost zero flow-rate. The high pressure is 
caused by the dynamic loss through the valves in the TPC flow cir¬ 
cuit that have been set to be almost closed. In these stalling 



(b) 


Fig. 5. Temperature recorded by (a) thermocouple TC-1, i.e. of the steam being generated in the HHX and flowing into the top of the displacement cylinder and (b) 
thermocouple TC-2 placed inside the CHX near the HHX-CHX connection, during normal steady-state operation. 
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conditions, the much reduced flow-rate leads gradually to a reduc¬ 
tion in the dynamic pressure loss, which allows the TPC to resume 
its pumping action, driving an increasing flow into the TPC flow 
circuit, increasing the pressure drop across the valves, and so on. 

Fig. 4 shows time plots equivalent to Fig. 3a, i.e. taken during 
normal steady-state operation, of the temperature of the solid 
walls of the F1HX (T H Hx-waii; in Fig. 4a) and of the CFIX (T C hx-w.iii; 
in Fig. 4b). Both plots belong to the same experimental run. The 
thermocouples from which the measurements of T H Hx-waii and 
IcHx-waii were obtained were positioned approximately 2/3 of the 
way into the HHX and CFIX, respectively. Similarly, Fig. 5 shows 
time plots of the working fluid temperature at the top of the 
HHX (r HH x-fiuid; Fig. 5a) and in the CHX (T CH x-nuid; Fig. 5b), at the 
same conditions. The location where these measurements were 
made are indicated in Fig. 1, which shows thermocouples TC-1 
(used for T HH x-fiuid) and TC-2 (used for T CH x-nuid). 

Fig. 4a reveals that the temperature of the HHX tubes T HH x-waii 
fluctuates between approximately 100 °C at the point in the cycle 
when HHX is flooded by liquid water (when y + is close to its peak 
value) and 550-600 °C when the HHX is empty of water (such that 
y~ is close to its peak value). Five stages are observed in the HHX 
temperature cycle. Starting from Point TV at t ~ 2 s in Fig. 4a: (i) 
a steep temperature rise from T HH x-waii ~ 100 °C to r HH x-waii 
~ 400 °C; (ii) a continuing temperature increase but with a smaller 
rate of change from T HH x-waii ~ 400 °C to the peak temperature 
FHHx-waii ~ 550-600 °C; (iii) a gradual temperature decrease from 
the peak temperature to T H Hx-waii ~ 400 °C; (iv) a steep tempera¬ 
ture drop down to T H Hx-waii ~ 100 °C; and (v) a short time spent 
at a constant temperature at the lowest end (T HH x-waii ~ 100 °C). 
The steep gradients in Stages (i and iv) are associated with the li¬ 
quid-vapour interface passing over the position of the thermocou¬ 
ple where the measurement is made. In Stage (i) the liquid water 
level passes this position in a downwards direction (d y + ldt< 0) as 
the water level falls from the HHX into the CHX. During this time, 
first the HHX and then the CHX are gradually emptied of water, 
which flows through the feedback line into the lower end of the 
power cylinder. The temperature of the HHX T HH x-waii rises sud¬ 
denly from 100 °C to 400-450 °C, as the solid surfaces become 
uncovered. Conversely, in Stage (iv) the water level passes the 
thermocouple position in an upwards direction (dy^/dt > 0), as 
the HHX is filled. The temperature of the HHX T HHX -waii drops 
abruptly from 400-450 °C to 100 °C as the solid surfaces are cov¬ 
ered by the water, and remains at this low temperature for as long 
as the water level is above the thermocouple location (this is Stage 
(v)). Stages (ii and iii) are recorded during the time when the water 
level is in the CHX and HHX respectively, and in both cases moving 
upwards (dy + /dt>0). Stage (iii) is related to the flow of steam 
through the HHX in an upwards direction. The steam is generated 
below the thermocouple position following the introduction of 
water from the CHX into the HHX, and then flows through HHX 
and the vapour connection into the top part of the power cylinder. 

The thermal stressing in the CHX walls is considerably smaller 
than that in the HHX. Fig. 4b shows a typical time trace of the 
CHX wall temperature, which indicates that this is of the order of 
30 °C, going from a low of 65 °C to a high of 95 °C in half the cycle 
period. It can be seen that, in the recorded conditions, no evapora¬ 
tion is possible in the CHX whose walls are retained at T c hx- 
wa ii < 100 °C by the flow of the cooling water on the tube side 
(see Fig. 1 and Section 2.1.1). 

The working fluid temperatures shown in Fig. 5 indicate that a 
small degree of superheating and subcooling takes place during the 
cycle. The saturation temperature of water/steam is not expected 
to change by more than a few degrees (typically 2-3 °C) due to 
the observed pressure fluctuations (~0.5 to 1 m). This is smaller 
than the observed 5-10 °C fluctuation in 7HHx-fiuid (~8 °C in 
Fig. 5a), which nonetheless is considered small. This is a useful 


finding as it confirms the conventional assumption that is made 
in the modelling of these systems [17-22], that the thermody¬ 
namic cycle occurs mainly within the saturation region without 
superheating or subcooling. 

3.2. Irregular operation 

Variability in the oscillatory behaviour of the TPC under consid¬ 
eration was observed during some experimental runs in the form 
of beating. The beating phenomenon was readily recognised both 
by visual observation and by an irregularity in the resulting pres¬ 
sure plots. An example of this is behaviour shown in Fig. 6. 

The observations of beating confirm earlier ones made in a sim¬ 
pler NIFTE prototype in Refs. [17,22], where the author linked the 
phenomenon to the large thermal capacitance of the heat exchan¬ 
ger blocks. Beating arises generally due to a slight mismatch be¬ 
tween two underlying simultaneous oscillatory signatures within 
the same signal (in our case, flow-rate or pressure). A possible rea¬ 
son for this observed behaviour in the present device could be the 
mismatching between the heat exchange and phase change pro¬ 
cesses taking part in the HHX and CHX, and the flow processes tak¬ 
ing part in the power cylinder and TPC flow circuit. 

Finally, in some experimental runs, especially with the short 
PTL, low amplitude oscillations were observed even when the 
valves connecting the pressure reservoirs were fully open, whence 
large amplitudes were expected. In some cases, the pump stopped 
pumping. It was not entirely clear what was causing this behav¬ 
iour; however, a possible reason for this could be a mismatch be¬ 
tween the HHX and CHX heat exchange capacities. 

It was noted in Refs. [17,21] that two loss mechanisms occur in 
this device due to the properties of the solid materials from which 
it is constructed, and that these mechanisms can lead to engine 
stall. The first mechanism, known as ‘entrance condensation’, is 
caused by the pressure variations within the engine and involves 
the periodic condensation and evaporation of the working fluid 
on the interior surfaces, which is accompanied by a periodic flow 
of heat into and out of the solid walls. As explained in Refs. 
[17,21], this leads to a direct loss in pressure during the discharge 
stroke as vapour condenses on the colder walls in the top part of 
the engine, and a parasitic pressure rise during the suction stroke 
as the condensate re-evaporates. The second loss mechanism, re¬ 
ferred to as the ‘shuttle loss’, occurs in the power cylinder over 
the height in which the liquid level moves upwards and down¬ 
wards, and involves the periodic storage of heat in the solid walls 
of the power cylinder during the high pressure discharge stroke, 
and the periodic removal of heat from the walls by the liquid dur¬ 
ing the low pressure suction stroke. The shuttle loss causes a direct 
pressure loss approximately in phase with the displacement of the 
liquid in the power cylinder. In addition, this loss mechanism was 
stated as being more pronounced than the entrance condensation 



Fig. 6. Pumping action with beating at two distinct frequencies. 
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(a) 


(b) 




(d) 


Fig. 7. Configuration A (small-capacity TPC) with (a) short = 2 m and (b) long = 5 m. Legends indicating the various settings for Q^,,! at which the U against AP results 
were obtained. Configuration B (large capacity HHX-CHX) with (c) short Lptl = 2 m and (d) long Lptl = 5 m. Legends indicating the various settings for Qboii at which the U 
against AP results were obtained. 


loss. Both loss mechanisms worsen with increasing load (i.e. pres¬ 
sure), which leads to increased temperature amplitudes in the va¬ 
pour containing regions of the engine. 

It is evident that these loss mechanisms occur predominantly at 
high loads and cannot be responsible for the stall observed at low 
loads in the present investigation. The finding of an additional loss 
mechanism in this device, and by extension similar oscillator en¬ 
gines, is an important development. Nevertheless, if the underlying 
reason for this behaviour is indeed linked to the design of the HHX 
and CHX as suggested above, we expect the solution to this prob¬ 
lem to be relatively straightforward. 

4. Results and discussion 

4.J. Pumping performance (U vs. AP) 

Regular operation reveals trends in which the mean rectified 
(unidirectional) flow-rate in the TPC flow circuit U decreases 
monotonically as the head AP that is applied across the TPC (i.e. 
between the HP and LP reservoir, with the use of the valves) is in¬ 
creased. This is also typical of common electrical circulators. As a 
consequence, the highest flow-rate U 0 (presented in Section 4.1.1) 
is achieved by the TPC when the head difference across the device 
is almost zero, and the maximum head achievable by the device 
AP max (presented in Section 4.1.2) is attained at stalling conditions 
when the flow-rate drops to ‘almost’ 2 zero. The U against AP perfor¬ 
mance characteristic of the TPC in the four different investigated 
configurations is summarised in Fig. 7. 

Fig. 7a shows results generated by the TPC with the small- 
capacity HHX-CHX combination and the short PTL (Lptl = 2 m). A 
similar pumping capability is observed at low AP between the 


2 Recall that AP is a dynamic (pressure) load across the TPC, which would also have 
fallen to zero had the flow through the CH and TPC flow circuits completely stalled. 


two employed heat settings Qboii = 10 and 15 kW. On the other 
hand, the stalling behaviour of the TPC (i.e. the operating point at 
near-zero flow-rate U and maximum head difference AP max ) dete¬ 
riorates at the higher heat setting (i.e. AP max reduces significantly). 
Fig. 7b concerns the same small-capacity TPC configuration as that 
in Fig. 7a, but with the PTL replaced by a longer length of hose 
(LpTL=5m). At low AP the flow-rate U is again similar at both 
Qboii = 10 and 15 kW heat settings, but once again, there is a signif¬ 
icant deterioration of the pumping capability at medium and high 
AP, including a greatly reduced AP max , as the boiler setting is in¬ 
creased from Qboii = 10 to 15 kW. 

Referring to Fig. 7c and d, the TPC configuration involving the 
large-capacity HHX-CHX combination shows a much more ro¬ 
bust high head performance, to the variations in the boiler heat 
setting, compared to the small-capacity TPC. With a short PTL of 
Lpt L = 2 m (Fig. 7c) there is an improved pumping performance at 
low AP due to the increase in Qboii, but a decrease in perfor¬ 
mance at higher AP, including a small reduction in AP max . Re¬ 
sults obtained with the longer LpT L =5m PTL (Fig. 7d) reveal a 
much more significant, and also non-monotonic, effect of Qboii 
on the pumping capacity at the lower to medium range of inves¬ 
tigated AP. Increasing Qboii from 10 to 12 kW results in a consid¬ 
erably reduced pumping performance for the same AP, but 
increasing Qb 0i i again from 12 to 15 kW leads to a slight 
improvement and recovery in U. On the other hand it is interest¬ 
ing to observe that, despite this earlier conclusion concerning 
the lower to medium APs, the maximum pumping flow-rate U 0 
is not noticeably affected by Qboii at the minimum (zero) head 
conditions. At the opposite end of the plots, i.e., high to stalling 
APs, the pumping flow-rates at the low and medium (Qboii = 10 
and 12 kW) heat settings are similar, but shows a small reduc¬ 
tion at higher Qboii = 15 kW, as was the case with the shorter 
PTL and also with the small-capacity TPC configurations dis¬ 
cussed in the previous paragraph. 
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Table 2 

Flow-rate U 0 in L/h at zero head AP = 0 m for the various pump configurations and 
Qboil- 


U 0 (L/h) in configuration 


Qboii (kW) 



HHX-CHX/TPC configuration 

Lptl (m) 

7 

10 

12 

15 

(A) Small-capacity 

2 

- 

650 

- 

650 


5 

- 

600 

- 

620 

(B) Large-capacity 

2 

550 

720 

720 

- 


5 

- 

850 

850 

850 


From an overall inspection of Fig. 7a-d it is possible to con¬ 
clude that the larger capacity (and size) heat exchangers of TPC 
Configuration B allow an improved performance (in particular, a 
higher flow-rate at a given AP) with respect the smaller Configu¬ 
ration A. This would have been expected, since a larger heat ex¬ 
changer area should normally lead to a greater heat exchanger, 
a greater volume of working fluid (water) involved in phase 
change, and thus, a greater ability to generate a positive displace¬ 
ment of the liquid piston into and out of the load. In addition, on 
comparing Fig. 7a-d, it is apparent that lengthening the PTL leads 
generally to an improved pumping capability for the TPC, irre¬ 
spective of the selected HHX-CHX combination. The only excep¬ 
tion to this (discussed in the next section) is found at zero 
loads with the small capacity HHX-CFIX combination. Certainly, 
a longer PTL is associated with an effective load, as seen by the 
TPC, that has both a higher flow resistance and a higher flow iner¬ 
tia. Given that the plots in Fig. 7 suggest that a higher flow resis¬ 
tance (and thus, higher dynamic head AP across the TPC) leads to 
lower flow-rate U, it may be concluded that the beneficial out¬ 
come of the longer PTL arises from the increased flow inertia in 
the load line, in agreement with the model predictions in Refs. 
[19,20,22], 


4.3.1. Maximum (zero head) flow-rate 

In Table 2 we compile the results for the maximum flow-rate U 0 
measured in zero head conditions for the various TPC configura¬ 
tions, extracted from Fig. 7. The results for U 0 are also shown in 
Fig. 8, as a function of the heat that is taken into the TPC at the 
HHX, Qhhx (see Fig. 1 for a definition of Qhhx. and also Table 1 
for the relation of Qhhx to Qboii). This heat would go towards evap¬ 
orating the liquid water within the TPC, which is expected to be the 
key process in determining the device’s pumping capability U 0 at 
zero load. Fig. 8 shows multiple results from a number of indepen¬ 
dent repeat runs (e.g., three separate runs with TPC Configuration B 
and Lpt L = 5 m at Qhhx ~ 3.6 kW), whereas Table 2 only contains an 
average value for U 0 at each investigated Qboii- 

As mentioned above, the best performance is associated with 
TPC Configuration B, i.e. the large-capacity TPC, and when 
employing the long PTL. It can be observed that in this configura¬ 
tion, U 0 is not affected by Qboii (or, Qhhx)- In fact, U 0 is not affected 
significantly by Qboii or Qhhx in any of the TPC configurations, ex¬ 
cept in the case of the large-capacity HFIX-CHX (Configuration B) 
with the short PTL, where it increases initially from 550 L/h at 
Qboii = 7 kW (Qhhx - 2.6 kW) to 720 L/h at 10 kW (Qhhx ~ 3.0 kW). 
It is noted however that, even in this case, a further increase in 
Qboii to 12 kW (Qhhx ~ 3.6 kW) produces no further strong in¬ 
crease in U 0 , although a very slight upwards trend may be dis¬ 
cerned in Fig. 8. Even with large variations in Qhhx (of a factor 
of ~2, going from about 2.5 kW to 5 kW), the pumping capability 
of a particular configuration at zero head U 0 does not show a 
noticeable increase, as would have been expected. Intuitively, 
higher heat inputs to the same TPC would be expected to yield 
increased quantities of steam from evaporation, increased dis¬ 
placements of the water in the power cylinder and consequently 


■S 


800 

600 

400 

200 

0 


S Configuration A.PTL=2m 

x X 

■ Configuration A,PTL=5m 


A Configuration B,PTL-2m 

▲ 

X Configuration B,PTL=5m 

• A A • 


▲ 


■ 


A j, A 

1 

2 3 4 


<W kW ] 


Fig. 8. Maximum flow-rate L! 0 at zero head for the various configurations as a 
function of Qhhx- 


increased flow rates. This interesting observation suggests that 
the HHXs in both TPC configurations possibly reach an optimum 
limit in their ability to transfer heat to (and thus to evaporate) 
the liquid water at Qhhx = 2.5-3.0 kW, after which only marginal 
increases in U 0 can come from further increases in Qhhx- The data 
in Fig. 8 indicate that for the large-capacity HHX this limit is pos¬ 
sibly at Qhhx ~ 3.0 kW. Moreover, it is expected that the small- 
capacity HHX would reach this limiting point at a lower heat 
flow-rate Qhhx. which may explain why this behaviour is not ob¬ 
served in the results from this TPC configuration, whereby U 0 is 
independent of Qboii and Qhhx within the experimental uncer¬ 
tainty and repeatability. 

The maximum pumping capacity U 0 of the small-capacity TCP 
configuration (in the approximate range 600-650 L/h over the 
range in Qhhx = 2.9-4.7 kW) lies in between 70% and 90% of that 
of the large-capacity TPC configuration (in the approximate range 
700-850 L/h over the range in Qhhx = 3.4-5.5 kW). The average 
discrepancy between the two pumps over the entire range of con¬ 
ditions is 75%, which is of the order of the difference in their avail¬ 
able area for heat exchange, which adds to the suggestion that the 
pumping performance of the TPC is strongly set by the HHX: the 
larger HHX is able to transfer larger quantities of heat from the boi¬ 
ler flame to evaporate water in the TPC, which results in increased 
flow rates at zero head. 

Furthermore, for a given TPC configuration the PTL also has an 
effect on the maximum pumping capacity U 0 over the investigated 
range of Qhhx- Lengthening the PTL from 2 to 5 m in the small- 
capacity TPC (Configuration A) results in a very marginal decrease 
in U 0 of about 6% or 30-50 L/h (from 650 to 600-620 L/h), whereas 
in the large-capacity TPC (Configuration B) the same lengthening of 
the PTL results in a more substantial increase in U 0 of about 20% or 
130-150 L/h, from 700-720 to 850 L/h. The decrease in U 0 in Con¬ 
figuration A can be explained to a large extent in terms of the ~10% 
decrease in the operating frequency of the TPC / 0 associated with 
the longer PTL (this is reported in Table 4 and discussed in Sec¬ 
tion 4.2.1). On the contrary, the decrease in U 0 in Configuration B 
cannot be explained in terms of / 0 , which in fact is unaffected by 
PTL (see Table 4), and can only be due to an increased capacity 
to evaporate the working fluid. It was stated above that a longer 
PTL is associated with a higher effective load flow inertia, which 
may be beneficial in the operation of the TPC [19,20,22], It is not 
clear for what reason this is not the case in the small-capacity 
TPC, especially given the fact that from Table 1 the heat into the 
HHX, Qhhx. at the same Qboii is approximately the same between 
the short and long PTLs. 

The observations above suggest that the heat input to the TPC, 
which should affect the maximum flow rate at zero head U 0 , is 
not determined by the heat setting Qboii, but by the design and 
choice of the HHX (and perhaps CHX, or other components), at 
least over the range of investigated conditions. Note that, ideally, 
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Table 3 

Stalling head AP max in m for the various configurations and Qb 0 n- 


APmax (m) in configuration 


Qboii (kW) 



HHX-CHX/TPC configuration 

Lptl (m) 

7 

10 

12 

15 

(A) Small-capacity 

2 

- 

5.0 

- 

3.2 


5 

- 

6.8 

- 

4.8 

(B) Large-capacity 

2 

7.0 

6.5 

5.9 

- 


5 

- 

8.4 

8.4 

7.5 


Table 4 

Frequency at zero head f 0 in Hz for the various configurations and Qboii. 


/o (Hz) in configuration 

Qboii (kW) 

HHX-CHX/TPC configuration 

Lrtl (m) 

10 12 15 

(A) Small-capacity 

2 

0.33 ± 0.01 

5 

0.29 ± 0.01 

(B) Large-capacity 

2 

0.24 ± 0.01 

5 


the heat input to the TPC would go towards evaporation and con¬ 
sequently displacement of fluid, especially at zero head. On com¬ 
paring the two TPC configurations with each other this effect 
should be considered alongside the fact that the frequency of the 
large-capacity TPC is lower by an average 25% compared to that 
achieved in the small-capacity TPC (again from Table 4). Hence, 
the 25% higher flow observed with the large-capacity TPC are ex¬ 
pected to arise as a consequence of increased volumetric steam 
generation ability of the order of ~60% to 70%. 

4.1.2. Maximum (stalling) head 

Table 3 and Fig. 9 summarise the maximum head capability of 
the TPC, AP milx , measured in stalling (negligible flow) conditions 
for the various configurations, again extracted from Fig. 7. From Ta¬ 
ble 3 and Fig. 9a it can be seen that, as with for maximum flow-rate 
U 0 , the best performance is associated with the large-capacity TPC 
and when employing the longer PTL. The improvement in AP milx 
when changing from the small capacity to the large-capacity TPC 
is between 25% and 55%, while the equivalent improvement when 
changing from the short to the long PTL is about 35-50% in Config¬ 
uration A and 30-40% in Configuration B. Specifically, these rela¬ 
tive improvements when changing the PTL correspond to 
noteworthy (and similar) increases in AP max by 1.5-2.5 mH 2 0 in 
both configurations, when increasing Lpj L from 2 to 5 m. 

Interestingly, AP max generally deteriorates (decreases) at higher 
Qboii, especially for the small-capacity TPC configuration where it 


drops by about 2 mH 2 0, e.g. from 6.8 m at lOkW to 4.8 m at 
15 kW with IpTL = 5m. In the large-capacity TPC configuration 
the decrease in AP max is still evident, but this is more limited, drop¬ 
ping by about 0.5-1.0 mH 2 0, e.g. from 8.4 m at 10 kW to 7.5 m at 
15 kW with LpT L = 5 m. In summary, the large-capacity TPC with a 
longer PTL and at the lower Qb 0i | setting has been demonstrated 
as being capable of operating with the maximum head differences 
across the circulator, which in the present experiments amounts to 
8.4 mH 2 0. 

Being a heat engine the TPC is expected to demonstrate a higher 
efficiency and a higher power output when the temperature differ¬ 
ence across it is increased. Thus, its maximum head capability is 
also expected to increase at higher temperature differences Ar in t. 
Here, AT int = T HH x-nuid - IcHx-nuid is the mean internal temperature 
difference experienced by the working fluid as it undergoes the 
thermodynamic cycle within the TPC, determined from tempera¬ 
ture recordings such as those shown in Fig. 5. Fig. 9b confirms that 
this is indeed the case for runs selected to be close to stalling at the 
highest heads, but not yet at fully stalled conditions. At full stall 
AT int reaches a maximum value of about 15-20 K, where it seems 
to saturate independently of AP max . The average gradient in the 
‘high pressure’ AP max per unit AT int in Fig. 9b is ~6 mH 2 0 per 7- 
8 K, suggesting that a good rule of thumb for engineering design re¬ 
quires an increase in Ar in t of the order of 1 K of a raising in AP max 
of 0.7 m. 

As a final comment, it was found that the maximum thermal 
efficiency attained by the TPC in this work was achieved at high 
head conditions, one setting prior to stall. The maximum value of 
the thermal efficiency amounted to about 2.5-3%, and was dis¬ 
played by the large-capacity TPC with the longer PTL and at the 
lower Qboii setting of 10 kW. This finding agrees well the conclu¬ 
sion above that this configuration and setting was also associated 
with the maximum head difference across the circulator. At zero 
head and maximum flow conditions the (useful) power output, 
and thus also, the thermal efficiency was zero. 

4.2. Dynamic behaviour (f and y) 

Thus far in the paper we have presented the measured pumping 
output of the TPC in its various configurations. In order to gain a 
better understating of the underlying processes that give rise to 
this pumping performance it is necessary to consider the operation 
of the TPC in more detail. The TPC is an oscillatory positive-dis¬ 
placement machine and hence its operational frequency / is a per¬ 
formance indicator of prime importance, as are the positive y + and 
negative displacement amplitudes with respect to the equilib¬ 
rium vertical position (defined in Fig. 1 ). These characteristics are 
the subject of the present section. 
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Fig. 9. Stalling pressure A P max for the various configurations as a function of (a) Qhhx and (b) AT int . 
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(a) (b) 


Fig. 10. Oscillation frequency/as a function of pressure difference AP for (a) Configuration A (small-capacity TPC) and (b) Configuration B (large-capacity TPC), with both the 
short Lptl = 2 m and long L PrL = 5 m. 




U [L/h] 

(a) 


Ap [mHiO] 

(b) 




(c) (d) 


Fig. 11. Oscillation amplitudes y* and y , and their ratio y + /y , as a function of (a and c) the mean flow-rate U and (b and d) the pressure difference AP, for Configuration A 
(small-capacity TPC) with the short Ipn. = 2 m. 


4.2.1. Oscillation frequency 

The frequency/of the TPC engine directly determines its output 
flow-rate, via the expression U = (y + - y~)Apf, where y + and y~ are 
the positive and negative displacement amplitudes of the liquid 
piston in the power cylinder (see Fig. 1) and A p is the cross-sec¬ 
tional area of the power cylinder. Fig 10 shows the variation in 
the measured frequency/of the TPC averaged over 10 experimen¬ 
tal realisations per run (fixed conditions) as a function of APfor the 
four different TPC configurations and at various Qb 0 n from 7 to 
15 kW. 

The oscillation frequencies / were found to be in the range 
0.30 ± 0.05 Hz in the small-capacity TPC configurations and 
0.20 ± 0.05 Hz in the large-capacity TPC configurations, decreasing 
slightly with increasing AP and almost independent of the setting 
of Qboii- The large-capacity TPC showed a very marginal decrease in 
/with AP, of the order of 0.01 Hz per 1 mH 2 0, whereas the large- 
capacity TPC had a more significant decrease of the order of 


0.04 Hz per 1 mH 2 0. Changing the length of the PTL from Lptl = 2 
to Lptl = 5 m resulted in a decrease in / by about 0.05 Hz when 
the small-capacity TPC configuration was tested. This finding is 
consistent with the explanation proposed earlier (in Section 4.1), 
that a longer PTL results in a greater liquid flow inertia in the load, 
which in turn would be expected to act so as to decrease the oscil¬ 
lation frequency of the TPC [19,20,22]. No consistent change in/ 
was observed when changing the PTL in the large-capacity TPC 
configuration. As with the flow-rate (L, the maximum frequency 
/o is observed at the lowest head AP («0 mH 2 0). 

The frequencies at zero head / 0 are summarised in Table 4. We 
note that the large-capacity TPC had a single / 0 of approximately 
0.24 Hz, independent of Qb 0 n and Lptl, while the small-capacity 
TPC/ 0 was slightly higher, ranging from 0.29 Hz (when the longer 
PTL was used) to 0.33 Hz (when the shorter PTL was used). For the 
same PTL, the increased liquid inertia and resistance that results 
from the longer length of piping in the large-capacity CHX would 
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Fig. 12. Oscillation amplitudes y + and y , and their ratio y + /y , as a function of (a and c) the mean flow-rate U and (b and d) the pressure difference AP, for Configuration A 
(small-capacity TPC) with the long Lptl = 5 m. 





(c) (d) 


Fig. 13. Oscillation amplitudes y* and y , and their ratio y'/y , as a function of (a and c) the mean flow-rate U and (b and d) the pressure difference AP, for Configuration B 
(large-capacity TPC) with the short L ptl = 2 m. 


also be expected to lead to reduced/[ 19,20,22 ], as indeed is the case. 
Lower frequencies are preferred for a given flow-rate capability, as 
this minimises the use of the check valves and cyclic thermal stres¬ 


ses in the HHX (and also, CHX and elsewhere). At least in this regard, 
the large-capacity TPC configuration with its lower/ 0 , but also gen¬ 
erally with its higher U for a given AP (see Fig. 7) is preferred. 
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U [L/h] 

(a) 





Ap [mH 2 0] 


(C) 


(d) 


Fig. 14. Oscillation amplitudes y* and y , and their ratio y*/y , as a function of (a and c) the mean flow-rate U and (b and d) the pressure difference AP, for Configuration B 
(large-capacity TPC) with the long Ipn. = 5 m. 


4.2.2. Oscillation amplitude (y) 

The output flow-rate of the TPC is also proportional to the dif¬ 
ference y* - between the positive (suction/inwards flow) y + 
and negative (discharge/outwards flow) and displacement 
amplitudes of the liquid piston in the power cylinder (see Fig. 1). 
Figs. 11 and 12 (for the small-capacity TPC) and Figs. 13 and 14 
(for the large-capacity TPC) show the absolute measured values 
of the amplitudes as a function of flow-rate U (a) and head rise 
across the TPC AP (b) for the various configurations. The same fig¬ 
ures also contain the relative ratio of the amplitudes y + /y” as a 
function of flow-rate U (c) and head rise AP (d). The common char¬ 
acteristic of these plots is that both y + and y~ increase (in absolute 
magnitude) in experiments in which the flow-rate U increases, 
whereas they decrease (in absolute magnitude) in experiments in 
which the head rise AP increases. This observation is consistent 
with the fact that the frequency / does not vary significantly for 
operation at different conditions, such that the output flow-rate 
is a strong function of the two displacement amplitudes. 

Proceeding further we observe that the ratio of suction to dis¬ 
charge amplitudes y + /y~ typically increases with U and decreases 
with head rise AP. The small-capacity TPC exhibits a ratio that in¬ 
creases approximately linearly from zero and reaches a maximum 
value of 0.6-0.8 at the highest U and lowest AP, with the upwards 
stroke being shorter than the downwards one. However, the large- 
capacity TPC has ratios exceeding unity, reaching values of 2.0 for 
Lpji = 2 m (Fig. 13) or 1.5 for L^l = 5 m (Fig. 14). It is possible that 
the greater pumping capacity of the large-capacity TPC is con¬ 
nected to these higher ratios, with the balance between the evap¬ 
oration and condensation ability of the HHX and CHX being the key 
determining factor in the relative negative and positive amplitudes 
exhibited by the TPC. Clearly, the sizing of the two heat exchangers 
is crucial in affecting this aspect of the device. The implication that 
results from this observation is that the improved capacity of the 
CHX (relative to the HHX) in the present study is what has allowed 
the improvements observed in Configuration B. 


5. Conclusions 

A thermally powered circulator based on the ‘Non-Inertive- 
Feedback Thermofluidic Engine’ (N1FTE), a two-phase thermoflui¬ 
dic oscillator/heat engine technology, was constructed and 
operated successfully. The circulator was employed as a 
replacement for a central heating hot water circulator coupled to 
a domestic gas-fired boiler. Two heat exchanger designs and two 
power transmission lines were tested; the latter being a length of 
tube filled with liquid (in this case, water) that connects the 
circulator hot water flow circuit (pumping load) to its power 
cylinder, within which an oscillating liquid piston displaces liquid 
into and out or the circulating hot water circuit. 

During regular operation the circulator achieved time-mean 
pumped flow-rates that decreased monotonically as the head ap¬ 
plied across it increased, similarly to common electrical circulators. 
The maximum flow-rate measured was 850 L/h, achieved with the 
high capacity combination of heat exchangers and a long power 
transmission line, whose increased inertia seems to offer perfor¬ 
mance improvements to the device. The maximum attainable 
(stalling) head measured was 8.4 mH 2 0, achieved with the same 
configuration. The maximum pumped flow-rates were attained at 
medium to high 10-15 kW boiler heat settings, and the maximum 
heads at medium 10-12 kW heat settings. Of this, the heat used by 
the circulator was about 3-5 kW. Further, it was found that an in¬ 
creased introduction of heat into the circulator led generally to 
greater pumping flow-rates, until a certain point, after which the 
flow-rates indicated saturation. Interestingly, the enhanced heat 
addition into the circulator also led eventually to a deterioration 
in the maximum pumping head capability. In addition, the internal 
temperature difference across the working fluid (water/steam) 
within the circulator showed a positive correlation with the max¬ 
imum head capability of the circulator, however, at the highest 
internal temperature differences of between 15 and 20 °C a limit 
in the maximum head performance was observed. 
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The pumped flow-rates achieved by the thermally powered cir¬ 
culator can be contrasted to the flow-rate requirement of approx¬ 
imately 500-700 L/h for a 10-15 kW central hearting boiler with 
60 °C supply and 40 °C hot water return temperature, respectively. 
This comparison implies that the prototype is capable of achieving 
a realistic pumping performance for this application, though it 
ought to be pointed out that at the higher heads that would expe¬ 
rienced in real central heating flow circuits the flow-rate output of 
the circulator deteriorates. 

The operational frequency of the oscillating circulator was not 
strongly affected by the conditions, and was set by the circulator 
configuration. It was found to be approximately 0.29-0.33 Hz for 
a circulator with short heat exchanger tubes, and 0.24 Hz for an 
alternate configuration with a longer length of exchanger tubes. 
In addition, a further increase in the length of the power transmis¬ 
sion line was shown to lead to lower operational frequencies, but 
slightly improved head and flow-rate performance. The latter con¬ 
clusion is in good agreement with predictions from simple dy¬ 
namic models of the circulator [18-22] that preceded the present 
study. Thus, it was found that the pumping capacity of the circula¬ 
tor was mostly determined by the oscillating liquid amplitudes in 
the power cylinder that defined the positive displacement ampli¬ 
tude of the liquid piston into and out of the hot water circuit. 
The best circulator configuration was associated with lower oper¬ 
ation frequencies and relatively large ratios of suction to discharge 
displacement. It is concluded that an improved circulator should 
result from an increased sizing of both the hot and cold heat 
exchangers to improve heat exchange capacity, but also from the 
careful balancing between the sizes of the two heat exchangers. 
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